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Abstract 

Recent results in the theory of multiphoton spectra for coherent radiation sources are overviewed, with the emphasis on channeling 
radiation. For the latter case, the importance of the order of resummation and averaging is illustrated. Limiting shapes of multipho¬ 
ton spectra at high intensity are discussed for different channeling regimes. In some spectral regions, there emerges an approximate 
correspondence between the radiative energy loss and the electron integrals of motion. 
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1. Introduction 

Ultrarelativistic electrons moving periodically in coherent ra¬ 
diation sources produce fairly monochromatic radiation. The 
intensity of this radiation is proportional to the length of the pe¬ 
riodic or conhning structure, and if it is sufficiently high, the 
electron is likely to emit several photons during its passage. 
The problem of multiple photon emission per electron plays 
s^nihcant role for crystalline radiators with thickness > l%Xo 
I 111 , Compton sources with optical laser pulse energy densities 
^ iQlsllj synchrotron radiation |3l, undulators with 

the number of periods > 10^ ||3l ■ At that, multiple photon emis¬ 
sions from different electrons can be discriminated by means of 
coincidences between detectors, except in cases when the elec¬ 
tron beam density is very high (e.g., in SASE FELs), and so 
only the problem of photon pileup from one electron remains. 

Physically, the multiple photon emission from an electron in 
a strong electromagnetic held may be viewed as a kind of elec¬ 
tromagnetic showering. It greatly simplihes if, as is typical for 
many coherent radiation sources, energies of emitted photons 
are much smaller than the electron energy: 

a> ^ E^. (1) 

Then, in the hrst approximation, it is permissible to count only 
photons stemming directly from the classical current generated 
by the initial electron, and neglect pairs produced by those 
photons downstream. 

For different experiments, though, the quantities of interest 
under conditions of multiple photon emission may differ. For 
instance, in synchrotrons one is concerned with lag Huctua- 
tions of an electron accelerated by the running wave, in FELs 
- electron bunching, in Compton sources - electron energy 
straggling before its last interaction with a photon yielding the 
observed gamma-quantum. For coherent radiation in oriented 
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crystals, one does not really care about the spent electron, but 
instead, measurements of the radiation are usually carried out 
with calorimeters, and are restricted to spectra integrated over 
(relativistically small) emission angles. Then, only the aggre¬ 
gate energy of all the nearly forward moving photons from one 
electron is measured. In the present article, we will discuss the 
latter case. 


2. Resummation technique and its alternatives 


As was mentioned above, the possibility of photon pileup 
in the calorimeter causes significant deviation of the multipho¬ 
ton probability (i.e., the radiative energy loss) spectrum dwfdcxj 
from the single-photon spectrum ^ predicted by 

classical electrodynamics, in spite of condition ([TJ. Nonethe¬ 
less, a multiphoton spectrum can be deduced from the corre¬ 
sponding single-photon one, as is described below. 

Technically, the situation is the simplest when one’s knowl¬ 
edge of initial conditions for the electron and the target is suffi¬ 
cient for completely specifying the electron’s trajectory. Then, 
presuming that photons exert negligible back action on the 
electron, their probability distribution obeys Poisson statistics 

iSIlIl 


dWn 

du\ . . . dtOn 


1 dw] dw] 
ni dci)\ dojn 


( 2 ) 


— r'” dto ^ 

where n is the number of emitted photons. Wo - e ■" ^ is 

the photon non-emission probability, and factor 1/n! accounts 
for photon equivalence. It is important to note, though, that 
the statistical independence of photon emission events does not 
imply that they occur sequentially in time. For instance, in a 
coherent state of the quantized electromagnetic field (a closest 


’Deviations from Poisson statistics may occur, e.g., in FELs, when a strong 
optical field seeded in the cavity stimulates emissions at the same frequency. 
But for gamma-ray sources, the transparency of any medium is too high to 
confine the photons and attain appreciable photon occupation numbers. 
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counterpart to a classical electromagnetic field) |0], there is a 
quantum superposition rather than mere statistical mixture of 
photon number states, which is indicative of their simultaneous 
generation. 

To derive from da the radiation spectrum measured by a 
calorimeter, one must integrate the fully differential spectrum 
for any definite number of emitted photons over their entire 
phase space, under the restriction of the sum of photon energies 
to a certain value at, and then sum over photon multiplicities: 


dw 
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OO ^ 
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Series (l3]l can further be resummecjl to all orders, by Laplace 
transforming it, whereafter it exponentiates, and then can be 
transformed back to the ai variable: 


(4) 

da> 2m Jc_,oo 

The integration contour in (|4]i is flexible in the complex s plane, 
but the integrand is not the only possible one, too: The upper 
limit of the integral in the exponent can be actually set equal to 
any value greater than a>, which can be expedient when deal¬ 
ing with incoherent radiation component 101 ■ For single-photon 
spectra with a relatively sharp upper end cjq <& Eg, it is usually 
sufficient to let the upper integration limit equal to cxJq. 

The practical benefit of the exponentiation is that it reduces 
the number of integrals to be handled. If one attempts instead 
to calculate termwise the probabilities of emission of a cer¬ 
tain number of photons via Eg. Q, computation of the mul¬ 
tiple integrals even numerically quickly becomes prohibitive. 
Sometimes, nonetheless, such multiple integrals are computed 
by Monte-Carlo - see, e.g.. Ref. 2 of IbIi. 

An approach alternative to (HJl is to view dwjdco as a solu¬ 
tion of an integro-differential equation with the single-photon 
spectrum playing the role of its kernel. Specifically, presuming 
that dwildtoi is proportional to the target thickness L, one can 
formally derive the standard transport equation 


d Q clw2 

—nCEe-w) = I dLOi — ^[n{Ee-a) + LO\)-TV{Eg-uj)] 
uL-/ uLi cioj\ 

(5) 

holding for the electron energy distribution function 


n(£, - w) = ^ + Wq{L)6{uj). 
dci) 


(6) 


Equation (|5]l might also serve as an alternative justification for 
the 1/n! factor in (|2]l, by virtue of the identity 
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^The term “resummation” was coined in the theory of radiative corrections 
at high-energy particle scattering (see, e.g., Ci). In view of the apparent tech¬ 
nical similarity, it seems to be proper also in the field of coherent radiation 
sources. 
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Figure 1: Simplified reconstruction of the single-photon spectrum from mul¬ 
tiphoton one, for radiation at volume reflection. Points, data of fT^ for a 2 
mm crystal bent to radius 11m. Blue band, fit of = Aci)e~^^ to exper¬ 
imental data, the adjusted parameters being A = 0.042 ± 0.003 GeV“^ and 
a = 0.053 ± 0.002 GeV“*. Red band, the reconstructed single-photon spectrum 
~ ^ ^ i-A/a^i_ Dashed green curve, prediction for the single¬ 

photon spectrum of coherent bremsstrahlung in a bent crystal B. 


should it not be the fact that each photon may actually be co¬ 
herently generated along a large length, simultaneously with 
formation of other photons. 

Another popular approach for computing multiphoton emis¬ 
sion probabilities during semi-classical electron passage is to 
divide the trajectory into several pieces, the interference of ra¬ 
diation from which is expected to be negligible, and define a 
probabilistic process of photon emission or non-emission on the 
junctions. It also permits taking into account incoherent scatter¬ 
ing during radiation, and other accompanying processes. That 
approach led to construction of Monte Carlo simulators widely 
used in modern practice iSlIIIllIll Till . At the present state of 
the art, such generators usually employ kernels of generation of 
only one photon, which is subsequently iterated - as well as in 
the resummation approach described above. 

Adopting an analogy with high-energy physics, one can 
envisage further development of simulators to include non- 
factorized evaluation of two-photon emission along with the 
appropriate account of virtual radiative corrections, simulation 
of a shower including coherent photon emission from 
pairs, etc. But for qualitative analysis, the resummation method 
within the region of its validity is advantageous by its simplic¬ 
ity. Below, we will confine ourselves to discussing a few exam¬ 
ples of its application. 

Before we proceed, it is worth pointing out that in case when 
there are no random parameters on which the spectrum signifi¬ 
cantly depends and which are to be averaged over, relation (HJl 
can be inverted to reconstruct the single-photon spectrum from 
the observed multiphoton one fl. In Eig. [T]we show an ex¬ 
ample where a single-photon spectrum of radiation at volume 
reflection was reconstructed (in a simplified fashion) from the 
experimental data. Without taking multiphoton effects into ac¬ 
count, the agreement between the experiment and the theory 
leaves much to be desired fl^ . but after the reconstruction it 
becomes acceptable. Needless to say, the reconstruction pro¬ 
cedure is expected to be efficient only if the radiation intensity 
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(photon multiplicity) is not too high. In what follows, we will 
discuss cases of high intensity and show that the correspond¬ 
ing spectra may have diverse shapes depending on the physical 
conditions. 


our framework, the anomaly must keep weak, i.e. its index 
a = L/Xq must remain small, or else the shower will not be 
dominated by soft photons. 


3. High-intensity limit and the departures from Gaussian- 
ity 


For pure coherent radiation, typically characterized by an 
edgy single-photon spectrum, multiple photon emission imi¬ 
tates the appearance of higher spectral harmonics. With the 
increase of the radiation intensity, such quasi-harmonics prolif¬ 
erate and coalesce into universal forms depending on just few 
parameters. For a definite and smooth electron trajectorjH, the 
(central) limiting form of the multiphoton spectrum is a Gaus¬ 
sian, as one would expect^ 


dwc 

da> 


1 



exp 


(oj - 

2wj 


(7) 


The corresponding regime may be physically categorized as a 
normal random walk. At large o) beyond the maximum, the 
spectrum decreases closer to an ordinary exponential rather 
than Gaussian law fl. That can be interpreted in terms of 
a maximal-step walk, at which the probability is a product of 
probabilities for each maximal step, the number of steps being 
~ w/wo. 

A correction to (|7]i breaking the scaling was obtained in . 
It engages the third spectral moment related to skewedness. Be¬ 
sides rendering the asymmetry to the spectrum, it slightly red- 
shifts its maximum. 

Another source of departure from the Gaussianity, even in the 
absence of dependence of the underlying single-photon spec¬ 
trum on external parameters, is the possible infinity of the spec¬ 
tral moments defining the Gaussian. For instance, in the case 
of coherent bremsstrahlung (at a highly over-barrier electron 
motion), the dependence on electron impact parameters is neg¬ 
ligible, but yet there is an incoherent radiation component de¬ 
scribed by Bethe-Heitler single-photon spectrum ^ for 

which the spectral moments diverge in the ultraviolet (or more 
precisely, are dominated by the upper endpoint oj Eg, around 
which our soft-photon resummation approach breaks down). 
That leads to a modification of the limiting spectrum shape by a 
power-law tail persisting at high u. The tail, however, does not 
override the whole spectrum, whose maximum is still mainly 
determined by the coherent part of the spectrum and its con¬ 
vergent first moments. Therefore, the diffusion regime may 
be regarded as weakly anomalous, and the result!^ distribu¬ 
tion as intermediate between Gaussian and Levy [SI]. Within 


^Such conditions can be rather well met in undulator radiation or Compton 
sources. 

'^In 0, moments ZiJf = fg dwicji and = fg for the 

coherent radiation component were designated as ZZTfj and Here, for 
simplicity we will not be considering the incoherent contribution at all, there¬ 
fore omitting subscript c. 


4. Channeling radiation: the impact of initial state averag¬ 
ing 


Our formulation so far held for radiation from a completely 
prescribed electromagnetic current. Next, it is important to ask 
what happens in case of channeling radiation, when the radi¬ 
ation intensity, and therewith the photon multiplicity, strongly 
depends on the electron/positron oscillation amplitude in the 
channel, in turn depending on the electron impact parameter 
and the angle of entrance to the crystal. Clearly, in the latter 
case the multiphoton spectrum must be duly averaged over the 
initial electron beam. 

An important point here is that the averaging must be carried 
out after the resummation, and this order of operations is signif¬ 
icant under the condition of high photon multiplicity. Clearly, 
in Eq. (|4|i, the average of an exponential does not reduce to an 
exponentiated average. 

Similarly, in Eq. (|5]l. 


r[(Eg — OJ + oji) 


d dw\ 
dL doJi 


^ — OJ + tui)) 


d dw\ 
dL doJi 


(9) 

invalidating the transport equation for the averaged distribution 
function with the averaged kernel. In fact, such equations were 
heuristically used in the literature both for overbarrier motion 
and for channeling conditions. At that, sometimes they in¬ 
volved adjustable parameters to be determined from fits to ob¬ 
servable spectra. However, since the fraction of particles near 
the bottom of the channel radiates weakly, but contributes to the 
total number of passed electrons strongly, estimates even of the 
mean photon multiplicity per electron based on kinetic equation 
with an averaged kernel can be deceptive. 

A related issue is that the knowledge of an averaged single¬ 
photon spectrum without detailing its dependence on the elec¬ 
tron initial conditions is insufficient for deducing the (averaged) 
multiphoton spectrum. Indeed, the averaged single-photon 
spectrum depends on a single variable oji, whereas 

the averaged multiphoton spectrum nontrivially depends on two 
variables: oj and the target thickness L, while the dependence 
on L is intertwined with the dependence on the electron initial 
conditions. It should be stressed yet that the dependence of the 
radiation on initial conditions must include both the incidence 
angle and the impact parameter. Whereas the dependence on 
the former is measurable in principle (cf., e.g., ilTT ). the de¬ 
pendence on the latter is not (there are no relativistic beams of 
angstrom transverse size), and thus inevitably must be modeled 
theoretically. 

Let us now discuss typical high-intensity channeling radia¬ 
tion spectrum shapes, for different types of confining wells, and 
for simplicity neglecting dechanneling. 
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4.1. Positron channeling 

The averaged multiphoton spectrum is most easily tractable 
for channeling of positrons, when the continuous potential 
of aligned axes or planes may be regarded as approximately 
parabolic in the entire channel (as is known to be the case for 
silicon crystals in orientation (110), (111), or (100)). Then, 
the particle trajectories in the channels are given by harmonic 
curves, and the single-photon radiation spectrum is explicitly 
expressible, too. It is proportional to the amplitude squared 
of the transverse motion, and therefore to the electron trans¬ 
verse energy E_i_, whereas the differential of the phase space 
over which the initial conditions are uniformly distributed can 
be written as dE^, with D the channeling well dimensionality: 

1 nE±max 

(...)^-^ dEf...^D (10) 

Here value 0-2 corresponds to capture to axial channeling 
from a beam with angular divergence > 0cEI H - 1 - to capture 
to planar channeling from a beam with angular divergence > Oc, 
or to axial channeling from a beam co-aligned with the axis and 
having divergence ^ 6c, and O - Ijl corresponds to planar 
channeling from a perfectly aligned beam with divergence 
Oc 

The analysis of the contour integral (|4]l averaged in this way 
is still rather involved, but at high intensity, one can view the 
spectrum just as a superposition of Gaussians. Inserting Q into 
Eq. (ITOt then yields 


dWc* O f' [" - ^(^rmax)]^ 

1 ) 

( 11 ) 

The latter integral can be expressed in terms of error functions, 
but we will restrict ourselves to elucidating its qualitative be¬ 
havior. 

At w < Err(£'j_max) - A, with A » -^wf(£’j.max), the maxi¬ 
mum of the integrand colcdI(E_i_!nnx) falls well within the 

integration interval, thereby giving the dominant contribution to 
the integral. Neglecting contributions from the endpoints leads 
to the result 


dwc+ 


D 


to 
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dto 
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( 12 ) 


It is independent of the variance Wj(£’j.jnax), and might be 
obtained by replacing partial Gaussians (including their pre¬ 
exponential factor) in (fTTb by the infinitesimally narrow dis¬ 
tribution 

6 [cO- ^(Eimax)] ■ 

That also demonstrates the existence of an approximate one-to- 
one correspondence between to and the transverse energy of the 
positron quantified by the ratio ^ = £±/£j.max- 


^In that case, there arises a practical question how to disentangle channel- 
ing radiation we are concerned with from that from overbarrier particles. One 
possibility is to use slightly bent crystals, in which channeled particles consid¬ 
erably deflect but do not significantly change their radiation spectrum, and then 
select radiation events in coincidence with deflected particles — see, e.g., 0 . 





(a) 


(o/aJiiE^max) 



(b) 


Figure 2: (a) High-intensity approximations HU for multiphoton channeling 

radiation from positrons, for (rt)D = «max ~ ^^(^±max)/‘^i(£^±max) = 10- 
Dashed curve, planar channeling radiation from monokinetic beam, D = 1/2. 
Solid curve, planar channeling radiation, £) = 1. Dot-dashed curve, axial chan¬ 
neling radiation, !D = 2. (b) High-intensity approximations for multiphoton 
channeling radiation from axially channeled electrons. Solid curve, Eq. USE 
for A„^ax)/'^J(7’max, Amax) = 10- Dashed curve, the same with a 

lower cutoff imin = '^min = 0-2. 


On the contrary, beyond the break, the integrand exponen¬ 
tially increases throughout the integration interval. Therefore, 
the main contribution to the integral comes from a vicinity of 
the endpoint f = 1, around which the ^-dependence of the ex¬ 
ponent can be linearized, whereas be put to unity, giving 


dWe+ D ■^2a»j(£’_Lmax) - 

~ ^ 2wj(A_Lniax) 

dto w>£3r(£±max)+A sjn ap- - wr^(£’±max) 


(13) 


Examples of behavior of function (fTTll are shown in Eig.|2a). 
The spectra strongly depend on the value of O, especially at 
low CO, in accord with asymptotics (fT2l i. That illustrates the 
profound effect of the averaging on multiphoton channeling ra¬ 
diation spectra. 


4.2. Axial channeling of electrons 

Compared to positrons, channeling of electrons, generally, 
is affected by dechanneling stronger, because electrons are at¬ 
tracted to atomic nuclei, and thus may enter regions of strong 
scattering. Eor axial orientation of the crystal, though, elec¬ 
trons can orbit around positively charged atomic strings, and 
thus keep on distance from nuclear vicinities. The approximate 
axial symmetry of a string alone, without introducing oversim¬ 
plified assumptions about the rj_ dependence, does not yet per¬ 
mit explicitly expressing electron trajectories and therefrom the 
radiation spectra. Nonetheless, qualitative analysis is feasible 
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for a potential modeled by a pure logarithmic function with a 
sharp cutoff at large radial distance 

V(r_L) = Vo In —- r^), < r,,. (14) 

ra 

The motion of electrons in potential (fT4l i is determined by 
two integrals of motion: transverse energy E_i_ and projection 
A of the angular momentum on the string direction. Estimates 
of first moments of the single-photon spectrum in a crystal of 


definite 

: thickness 

give: 
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where T is the period of the radial motion, exponentially de¬ 
pending on Ei_IVo, and relatively weakly dependent on A. Both 
spectral moments blow up as T ^ 0 or A —> 0, i.e., when the 
electron passes close to the singularity of the continuous poten¬ 
tial. 

The element of the uniformly populated phase space can be 
expressed in terms of T, A, as well: 

d^rj_d^p± oc TdEj_dA ^ VodTdA, (17) 


and the region accessible to channeling (subject to condition 
< oil is 


0 < d = 


|A| ^ T 

- < T = - 


< 1 , 


where r^ax ~ yfEJV^, A^ax = fa ^EeVole. Ultimately, the 
averaged multiphoton spectrum assumes the form 


^/Wg-ax 

dco 


2 



) 



[cU/tr U>i(^niax)] 
2/lw^(2niax) 


, (18) 


with 2max = TmaxAmax- Again, we restrict ourselves to qualita¬ 
tive analysis of integral (fTSl l. 

At w < wr(2max) - A, A » ■\fo^{Qmax), the maximum of 
the integrand at At = coUQ^nux)/<^ is beyond the integration 
domain. Therefore, the main contribution comes from the end¬ 
point A - T - 1, and generally, the integral is exponentially 
small: 


dWg-^x 2 oj^ (2max) 

doj -<t5r(G„„)-A V ;r - a;2] [ZuKgmax) + w] 


X exp 


[cO tUi(^niax)] 
2w^(2max) 


(19) 


^Here we refrain ourselves from consideration of channeled particles with 
> 0, which exist for axial channeling. The contribution from such particles 
is relatively small and will be analyzed elsewhere. 


On the contrary, at w > arr(2niax) + A, the maximum of the 
integrand falls well within the integration domain. Neglecting 
the endpoint contributions, or simply inserting 

b [W - wr(2max)/dT] (20) 


in place of the partial Gaussians, the asymptotics of the multi¬ 
photon spectrum derives as 


dWg-six (^liQmax) , dl 

- - -^-In^-. 

da> oj>oJi(Q^.j^)+A (x)^ (^l(Qmax) 


( 21 ) 


The corresponding mean energy loss dojoj '^'^'’^" diverges 
in the ultraviolet as a double logarithm, in accordance with 
the double logarithmic divergence of the single-photon average 
mean energy loss, cf. Eq. 

The delta function (l20l i implies that a correspondence be¬ 
tween a> and the integrals of motion holds for electron axial 
channeling radiation, as well. But in contrast to radiation from 
channeled positrons, the correspondence here takes place at the 
outer slope of the spectrum [see Eig.|2jb)]. Yet, since there are 
two integrals of motion, the correspondence of co with them is 
no longer one-to-one, unless one refers to the special integral of 
motion Q - TAoc 


Experimentally measured multiphoton axial channeling ra¬ 
diation spectra llTl IlST decrease considerably faster than by 
power law (|2TJ. That must be attributed to the crudeness of 
our model (ITtIi and the physical finiteness of (uTT) due to the 
smearing of the string potential at low rx_ through thermal and 
quantum vibrations. The simplest way to take all that into ac¬ 
count is to introduce a lower cutoff dmin = Tmin into integrals 
(HI, which leads to an exponentially decreasing spectrum [see 
Eig.|2b), dashed curve]. The beginning of the exponential sup¬ 
pression region depends on the value of the cutoff parameter 
and therethrough on the crystal temperature. 


5. Limitations and challenges 

The intrinsic simplicity of the soft photon resummation pro¬ 
cedure makes it a useful tool in various physical situations, pre¬ 
dicting generation of Imher quasi-harmonics at moderate inten¬ 
sity (as discussed in ISD) and the existence of a set of limiting 
distributions at high intensity, depending on the shape of the 
confining well, yet yielding in certain spectral regions approx¬ 
imate correspondence between the radiative loss u) and one of 
the integrals of motion of the channeled particle. 

At the same time, the resummation method has its limitations 
and is unable to cover, e.g., the case when a considerable part 
of the spectrum concentrates near the kinematic edge to - Eg 
10. Also, even in the domain of applicability of the soft pho¬ 
ton approximation ([T]), one may wish to evaluate higher order 
corrections in w/Eg, e.g., in order to account for degradation of 
Eg stipulating additional broadening and redshift of the multi¬ 
photon spectra yOfl, or to describe the effect of electron beam 
bunching in EELs fa. 

Eor the specific problem of channeling radiation, there arise 
additional issues. Eirst is the need for knowledge of not only the 
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averaged single-photon spectrum, which might be measured at 
a sufficiently small crystal thickness, but also of the dependence 
of the radiation spectrum on electron entrance parameters. That 
can be illustrated, e.g., by the variance of the multiphoton spec¬ 
trum about the mean value, for positron channeling in a har¬ 
monic well discussed in Sec. 14. II 



( 22 ) 


Due to the second term, it depends on the target thickness non- 
linearly, yet involving a dependence on the value of 2), which 
in turn depends on the shape of the well. Also, the existence of 
a correspondence in some spectral regions between u) and the 
particle integrals of motion implies that not all the values of in¬ 
tegrals of motion contribute at a given co, in contrast to the case 
of the averaged single-photon spectrum. 

Further complications emerge because in the presence of 
dechanneling, relative contributions of different integrals of 
motion can attenuate differently. The corresponding factors 
thus must be introduced into integrals (fTTb . (fTSl l. Thirdly, the 
dechanneled particles contribute to the radiation spectrum, as 
well, although in a somewhat harder spectral region. That de¬ 
mands inclusion of their contribution to the multiphoton spec¬ 
trum, too. On the other hand, studying multiphoton radiation 
spectra under such conditions may offer deeper insight into the 
complex dynamics of channeling. 
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